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1. Introduction 

Symmetry-breaking phase transitions are ubiquitous in condensed matter 
systems and in quantum field theories. There is also good reason to believe 
that they feature in the very early history of the Universe. At many such 
transitions topological defects of one kind or another are formed. Because of 
their inherent stability, they can have important effects on the subsequent 
behaviour of the system. 

In the first of these lectures I shall review a number of examples of 
spontaneous symmetry breaking, many of which will be discussed in more 
detail by other lecturers, and discuss their general features. The second 
lecture will be mainly devoted to the conditions under which topological 
defects can appear and their classification in terms of homotopy groups 
of the underlying vacuum manifold. In my final lecture, I will discuss the 
'cosmology in the laboratory' experiments which have been done to try 
to test some of the ideas thrown up by discussions of defect formation in 
the early Universe by looking at analogous processes in condensed-matter 
systems. 

2. Spontaneous symmetry breaking 

Often a system has symmetries that are not shared by its ground state or 
vacuum state. This is the phenomenon of spontaneous symmetry breaking. 
It always implies a degeneracy of the ground state. In this lecture, I want 
to discuss a number of simple examples of this phenomenon. 

The key signature of spontaneous symmetry breaking is the existence 
of some operator cf>, the order parameter (for example the magnetization 
M in a ferromagnet) whose ground-state expectation value is not invari- 
ant. Typically, the equilibrium state at high temperature is invariant - 
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for example, above the Curie temperature T c the expectation value of M 
vanishes. Symmetry breaking appears as the system is cooled below T c . It 
spontaneously acquires magnetization in some random direction. 

2.1. FERROMAGNET 

As a first example, let us think about a Heisenberg ferromagnet [1], a system 
of spins S r at lattice sites r, with Hamiltonian 

6 = ~\ E E J d r " r 'DSr Sr' - H £ Sr, (1) 
r r ' r 

where J, assumed positive, is non-zero only for neighbouring spins, and H 
represents a possible external magnetic field. For non-zero H, the ground 
state |0) of the system has all spins aligned in the direction of H: 

(0|S r |0) = lh-a, h = ^. (2) 

If we take the infinite-volume limit and then let H —> 0, this remains true. 

Clearly, with H = 0, (1) is invariant under all rotations. The symmetry 
group is 50(3) or, rather, its two-fold covering group SU(2). However, 
(0|S r |0) is evidently not invariant. Applying the rotation operators to |0) 
yields an infinitely degenerate set of ground states | Oh) , labelled by the 
directions of the unit vector h. Note that the scalar product of two distinct 
ground states tends to zero in the infinite- volume limit: 

|(0 h |0 h ')| = cos JV |-»0 as N^oo, (3) 

where A?" is the number of lattice sites and 9 is the angle between h and 
h'. In fact, these two states | Oh) and |0h') belong to separate, mutually 
orthogonal Hilbert spaces, with unitarily inequivalent representations of 
the commutation relations. No operator constructed from a finite set of 
spins has a nonvanishing matrix element between |0h) and |0h')- This is 
another characteristic of spontaneously broken symmetry. 

Note that physically, unless we introduce a symmetry-breaking inter- 
action with an external magnetic field, the different ground states are in- 
distinguishable. Indeed it is possible (though not particularly helpful) to 
define an invariant ground state, a uniform superposition of all |0h). 
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2.2. FREE BOSE GAS 

Another example of spontaneous symmetry breaking is Bose-Einstein con- 
densation. Let us first consider a free Bose gas [2], with Hamiltonian 

H = ^ e ka k «k, (4) 

k 

where a k and a k are the destruction and creation operators for a particle 
of momentum Kk, satisfying the commutation relations 

[a k ,a k ,] = <5 k;k /, (5) 

and 

£k = ^T- (6) 

The relevant symmetry here is the phase symmetry 

a k i ^ e^a k , a k i ^ e _M? a k , (7) 
corresponding to the existence of a conserved particle number, 

iV = E%«k. (8) 

k 

In the grand canonical ensemble, the mean occupation number of mode 
k is given by the Bose-Einstein distribution 

(«L«k) = iVk = e ^} ll) _ r (9) 

where (3 = l/k&T is the inverse temperature and \i the chemical potential. 
If the mean number density n = (N) /CI (where fl is the volume) is specified, 
then the value of /i for any large temperature T is given by 

_ 1 V 1 f d 3 k 1 

n Y e^(^-M) - 1 ~ J (2tt) 3 e /3(^"^) _ i • [W> 

As T falls, jJL increases. Clearly, to avoid a divergence, fi can never become 
positive, but eventually it approaches zero at some critical temperature 
T c . At that point, the occupation number (9) diverges at k = 0, so the 
integral approximation in (10) is no longer adequate. We have to separate 
the k = mode from the rest. For all the other values of k, we can still 
make the integral approximation, and set \i = 0. Thus we get 



(11) 
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For T <T C , (i becomes extremely small, of order —IcsT/nO,. Thus, a macro- 
scopically significant fraction of all the particles is to be found in the single 
zero-momentum mode. 

Note that from (11) we can calculate the critical temperature T c at 
which no first becomes large. It is given by 



1 (2mk B T c \ 3 ' 2 f^x^dx 
n = T^Zyi *2 L ( 12 ) 



which yields 



(2vr) 2 V h 2 J Jo e* - 1 ' 



2vr h 2 n 2 ^ h 2 n 2 / 3 
Tc= C(|) 2 / 3 mkB =3 - 31 ^7' (13) 



2.3. DEGENERATE GROUND STATE 

To see how a degenerate ground state emerges in this case, it is more 
convenient to work in terms of the scalar field 

<M r ) = ^E«ke* k ' r . (14) 

k 

Then we can write 

r % 2 

H= d 3 r — V0 f (r) • V0(r). (15) 
Jn 2m 

The phase symmetry (7) can now be expressed as 

0(r) i ^ e^^r), ${t) ^ e~^^{r). (16) 

Now let us introduce an explicit symmetry-breaking source term playing 
the same role as an external magnetic field for the ferromagnet. In this case, 
we take 

H ± = - [ d 3 r[j*j>(r) + 0(r)} = -n 1 ' 2 (j*&o+ A), (17) 
Jn 

where j is a complex number, and consider the limit j ' — ► only after 
letting O — > co. Here Hi may be thought of as representing the possibility 
of particle exchange with the environment. We then find that the density 
operator po for the zero-momentum mode is the same as with j = but 
with a shifted field: 

po = (1 - e^e^o j ( 18 ) 

where 

b = a + J -. (19) 
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To find [X we must again use (10) and (11) but appropriately modified, 
namely 

' " (20 ) 



" e-^-1 /i 2 

Clearly, as f2 — > oo, the first term on the right becomes negligible, and so 

(21) 

(We assume for simplicity that \j\ <C y/no/fi, so that |U remains small 
enough to be negligible for the k / modes.) Then we find that the 
ground state for the zero-momentum mode is a coherent state \<p), such 
that bo\<p) = 0, in which the expectation value of is non-zero, with a 
phase constrained by the phase of j: 

(fa))=4> = rp» = r& /2 -L. (22) 

The coherent states may be labelled by the value of 4>, and are given 
explicitly by 

\<f>) =exp[f} 1 / 2 (,Mo-<fao)]|0). (23) 

Once again, we have a set of degenerate ground states labelled by a phase 
angle None of these states has a definite particle number. However, as 
— > oo the uncertainty in no tends to zero. Moreover, as before, the scalar 
product between states with different values of ■& tends to zero; one finds 

fae^e**) = exp[-lV(l-e- i{,? ~ x) )] ^0 as fi oo. (24) 

Here too, these states belong to distinct orthogonal Hilbert spaces, car- 
rying unitarily inequivalent representations of the canonical commutation 
relations. 

It is important to note that degenerate ground states and symmetry 
breaking occur only in the infinite-volume limit. In a finite volume, there 
is always a unique ground state, a uniform superposition of 'ground states' 
with different phases. As in the case of the ferromagnet, the limits of infinite 
volume and zero external field do not commute. If we keep j finite and let 
$7 — > oo, we arrive at a state with (4>(r)) = r/e 1 ^. On the other hand if we 
let j — > first, the particle number remains definite, and we end up with a 
superposition of all the degenerate ground states. Physically (unless j / 0) 
the states are indistinguishable. 
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A Bose gas with interaction is described by the Hamiltonian 

r h 2 
H = / d 3 r — V^(r) • V0(r) 
Jn 2m 

+1 f d 3 r f dV V{\r' -r\)ft(r')ft(r)j>{r)0(r'). (25) 
2 Jn Jn 

which is clearly still invariant under the phase transformations (16). So 
long as the interaction is weak, the qualitative picture is largely unchanged. 
Bose-Einstein condensation has been observed in alkali metal vapours, such 
as rubidium [3]. 



2.4. LIQUID HELIUM-4 

The normal-to-superfluid 'lambda transition' at a temperature of about 2 K 
in liquid 4 He is another example of a symmetry-breaking phase transitions 
in a bosonic system [4]. But unlike the case of alkali metal vapours, this is a 
system with strong interparticle interactions that substantially change the 
picture. Nevertheless, there are close similarities. The broken symmetry is 
still the phase symmetry associated with conservation of particle number; 
the transition may still be described in terms of a scalar field <p which 
acquires a non-zero expectation value below the transition. 

The most obvious difference caused by the interatomic interactions is 
in the nature of the excitations above the ground state. In 4 He at low 
temperature, these are collective excitations, phonons, whose dispersion 
relation is linear near k = 0: 



e k = c s h\k\, (26) 

where c s is the sound speed. At larger values of |k| the graph curves down- 
wards to a minimum; excitations near the minimum are called rotons. 

A useful description of 4 He is provided by the Ginzburg-Landau model 
[4], which may also be applied to other systems. The starting point is to 
consider the free energy, F say, as a function of the order parameter <p. 
At least in principle, F[<fi] may be calculated in the usual way from the 
partition function Z, by restricting the sum over states to those with a 
given expectation value of the order parameter field, (0(r)) = (f>(r). In the 
neighbourhood of the transition temperature, F can be expanded in powers 
of 0: 



F W\ = J^r^-V^^-V^r) 

+ / d\ [a(T)|0(r)| 2 + i/3(T)|<Xr)| 4 + ■ ■ -1 , (27) 
Jn 1 J 
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where the higher terms are usually unimportant, at least for a qualitative 
description. The coefficient (3(T) is always positive and may usually be 
taken to be constant. At high temperature, a(T) is also positive so that 
the minimum of the free energy occurs at <j) = 0. At low temperature, 
however, a becomes negative, and the minimum occurs at 



-a(T)\ 



1/2 



P(T) 



) • (28) 



As usual, the phase of (f) is arbitrary: we have a degenerate equilibrium 
state. 

Note that the critical temperature T c is the temperature at which a(T) 
= 0. Close to that point, we may take 

a(T)wai-(T-r c ). (29) 

A good qualitative picture of the behaviour of 4 He is given by the two- 
fluid model, normal plus superfluid. The scalar field (j) = {<fi) describes the 
superfluid component, defining both the superfluid density and velocity: 

n s = \4>\ 2 , v s = ^-W, where <f> = \<f>\e M . (30) 

The normal component corresponds to single-particle (or, rather, single- 
quasiparticle) excitations above the ground state. 



2.5. SUPERCONDUCTORS 

The electrons in a solid constitute a Fermi gas rather than a Bose gas. It 
is not single electrons that condense but bound pairs of electrons, Cooper 
pairs [4]. There is an effective attractive force between electrons near the 
Fermi surface k 2 = fcf. At least for conventional superconductors, this force 
preferentially binds pairs with equal and opposite momenta and spins. 
Below the critical temperature, we find that in the ground state 

(a kT a_ ki ) = F(k) / 0, for k = |k| ~ k F . (31) 

The order parameter <j> in this case can be taken to be an integral over such 
products of pairs of destruction operators, multiplied by the internal wave 
function of a Cooper pair. 

There is an important difference between this and the examples dis- 
cussed previously. The symmetry here is again the phase symmetry (16), 
but it is now a local, gauge symmetry: $ is allowed to be a function of space 
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and time, r). This is possible because of the coupling to the electro- 
magnetic field A fM (x) which transforms as 

^)=^)-^(4 (32) 

The factor of 2e in the denominator appears because this is the charge of 
a Cooper pair. It ensures that the covariant derivative 

Dj = + 2ijU^, (33) 

transforms in the same way as <j) itself. The Ginzburg-Landau model may 
be used for superconductors too, provided that the derivatives V(f> in (27) 
are replaced by covariant derivatives D</>. 

Symmetry breaking in gauge theories is a somewhat problematic con- 
cept. Indeed Elitzur's theorem [5, 6] says that spontaneous breaking of a 
local, gauge symmetry is impossible! — which might be thought to imply 
that what I have just told you is nonsense. More specifically, it says that, 
while for a global symmetry taking the infinite-volume limit and then let- 
ting j — > may yield a state with (0) / 0, in a gauge theory we always 
have 

lim lim ((f>) = 0. (34) 
j^o n— »oo 

But one must be careful not to misinterpret this (entirely correct) theorem. 
It applies only in an explictly gauge-invariant formalism. If, as is often done, 
we add a gauge-fixing term that explictly breaks the local symmetry (e.g., 
by imposing the Coulomb gauge condition V • A = 0) then the remaining 
global symmetry can be broken spontaneously. We certainly can define and 
use gauge-non-invariant states with (cp) / 0, though there must always be 
an alternative (but often inconvenient) gauge-invariant description. 

A model widely used as an exemplar of symmetry breaking in particle 
physics is the Abelian Higgs model, the relativistic version of the Ginzburg- 
Landau model. It is described by the action integral 

1 = Jd A x [-\F^ + D^*D^ - ±A(0"V - r] 2 ) 2 ] , (35) 

with 

D„<f> = + ieA^, = d^A u - d l/ A /Ji . (36) 

Here the coupling constant A plays the role of (3 and — Xrj 2 that of a. 
(Note that here, and in general when dealing with relativistic models, I set 
c = h = 1.) Symmetry breaking in this model is very similar to that in a 
superconductor. 
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2.6. LIQUID CRYSTALS 

A very different example is provided by the isotropic to nematic transition 
in a liquid crystal [7]. A nematic liquid crystal is typically composed of 
rod-shaped molecules that like to line up parallel to one another. There is 
no long-range translational order: the molecules are free to flow past one 
another. But there is long-range orientational order. At any point in the 
liquid there is a preferred direction, characterized by a unit vector n, the 
director. Note that n and — n are completely equivalent. 

The symmetry group here is the rotation group SO (3). Above the transi- 
tion temperature T c , the system is completely isotropic, with the molecules 
randomly oriented, but below it, the rotational symmetry is broken. 

A convenient choice for the order parameter in this case is the average 
mass quadrupole tensor Q of the molecules in a small region. When the 
directions of the molecules are isotropically distributed, Q = 0. But if they 
are aligned in the direction of n, it has the form 

Q = Q(3nn-l). (37) 

In particular, if n is in the z direction, then 

f-Q \ 
Q = -Q . (38) 
V 2Q) 



2.7. GENERIC CASE 

Let us now examine the generic situation. (For more detail, see for example 
[8].) Suppose the system has a symmetry group G. In other words, the 
Hamiltonian H is invariant under every operation g £ G: 

U L (g)HU (g) = H for all g £ G, (39) 

where U(g) is the unitary operator representing the operation g on the 
Hilbert space. 

However, we assume also that there is an operator 4> with a nonvanish- 
ing ground-state expectation value which transforms non-trivially under G. 
Specifically, we consider a multiplet of operators <p = (<f>i)i=i... n transform- 
ing according to some n-dimensional representation D of G: 

U- 1 (g)i i U(g)=J2D ij (g)j> j , (40) 

3 



or more concisely 

U-\g)4>U(g) = D(g)l 



(41) 
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We suppose that the expectation value in the ground state |0) 



(O|0|O) 



00, 



(42) 



say, is not invariant: 



(0\U-Hg)4U(g)\0) 



D(g)<h ± 0o 



(43) 



for some g € G. Obviously this implies that the ground state |0) is not 
invariant: 



the ground state is degenerate. 

In general, not all elements of G lead to distinct ground states. There 
may be some subgroup H of elements such that 



The distinct degenerate ground states correspond to the distinct values of 
<fi = D(g)<pQ. Hence they are in one-to-one correspondence with the left 
cosets of H in G (sets of elements of the form gH). These cosets are the 
elements of the quotient space 



This space may be regarded as the vacuum manifold or manifold of degen- 
erate ground states. 

For example, for a Heisenberg ferromagnet, G = SU(2), and H = U(l), 
the subgroup of rotations about the direction of the magnetization vector. 
Here M = SU(2)/U(l) = S 2 , a two-sphere. For a Bose gas, G = U(l), and 
H comprises the identity element only, H = 1 = {1}. Thus M = S 1 , the 
circle. 

A nematic is a slightly less trivial example. Here G = SO(3); however, H 
is not merely the subgroup SO (2) = of rotations about n. Rather, H is 
the infinite dihedral group, H = D^, which includes also rotations through 
7r about axes perpendicular to n. Correspondingly M is not the two-sphere 
but the real projective space RP 2 , obtained from S 2 by identification of 
opposite points. 

2.8. HELIUM-3 

Finally let me turn to the particularly interesting, and relatively compli- 
cated, case of 3 He [9, 10]. This lighter isotope also exhibits a phase transi- 
tion, though at a much lower temperature than 4 He, between 2 and 3 mK. 




D(h)<f>o = 0o for all h€H. 



(45) 



M = G/H. 



(46) 



11 



It is of course a Fermi liquid. So the mechanism of superfluidity is very 
different, similar to that of superconductivity. In this case we have Cooper 
pairs not of electrons but of 3 He atoms. The order parameter can again be 
constructed from pairs of destruction operators. 

There is however an important difference. In the original BCS model, 
the pairs were bound in an isotropic *S state; as indicated by the form of 
(31). But for a pair of 3 He atoms close to the Fermi surface it turns out 
that the most attractive state is the 3 P. The pairs have both unit orbital 
and unit spin angular momenta: L = S = 1. We need to consider a more 
general form of order parameter, related to the quantity 

F ab (k) = {a ka a- kb ), a, b =|, | . (47) 

The fact that S = 1 tells us that F should be symmetric in the spin indices 
a and b, so it can be expanded in terms of the three independent symmetric 
2x2 matrices, a^ian where <jj are the Pauli matrices. The fact that L = 1 
means that F should be proportional to k times a function of k = |k| only. 
Thus we can write 

F ab (k) =F(k)A lJ (a i ia 2 ) ab k j , (48) 

where the two-index tensor A may be normalized by tr(A^ J 4) = 1. 

The order parameter is essentially A times a scalar factor representing 
the density of Cooper pairs. Since it is now a 3 x 3 complex matrix rather 
than a scalar, the possible patterns of symmetry breaking are much more 
complex. There are in fact two distinct superfluid phases, 3 He-vl and 3 He- 
B, which are stable in different regions of the phase diagram; the A phase is 
stable only at high pressure and at temperatures not far below the critical 
temperature. In the presence of a magnetic field, there is a third stable 
phase, the Al phase. 

The system exhibits a much larger symmetry than 4 He. To a good ap- 
proximation, it is symmetric under independent orbital and spin rotations, 
as well as under the phase rotations as before. Thus the symmetry group 
is 

G = SO(3) s x SO(3) L x 17(1) (49) 

(There is also a weak spin-orbit coupling, whose effects I will discuss a little 
later.) 

In the A phase, the order parameter takes the form 

Aij = —=di{m,j + inj), (50) 
y2 

where d, m and n are unit vectors, with m _L n. The vector d defines an 
axis along which the component of S vanishes. If we define 1 = m A n, 
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then 1 is an axis along which the component of L is +1. In this case, the 
subgroup H that leaves A invariant comprises spin rotations about the 
direction of d, orbital rotations about 1 combined with compensating phase 
transformations, and, finally, the discrete transformation that reverses the 
signs of all three vectors. Hence 

H A = 17(1) x U(l) x Z 2 . (51) 

Correspondingly, the vacuum manifold is 

M A = G/H A = S 2 x 50(3)/Z 2 . (52) 

Here, the elements of S 2 label the direction of d, while SO (3) describes the 
orientation of the orthonormal triad (1, m, n). The Z 2 factor represents the 
identification (d, m, n) = (— d, — m, — n). 

The B phase, by contrast, is characterized by an order parameter of the 
form 

Aij = Rije , (53) 

where R £ 50(3) is a real, orthogonal matrix. In this case, the only ele- 
ments of H are combined orbital and spin rotations, so 

H B = SO(3) and M B = G/H B = SO(3) x S 1 . (54) 

As I mentioned earlier, there is actually a weak spin-orbit coupling term 
in the Hamiltonian, which is only noticeable at long range, and which re- 
duces the symmetry to 

G' = SO{3) 3 x 17(1) with J = L + S. (55) 

Note that going from G to G' is not strictly speaking a case of spontaneous 
symmetry breaking. There are similarities: at short range, the symmetry 
appears to be the larger group G; when we go to long range (or low energy), 
we see that the symmetry group is actually G' . However, the true symmetry 
is always G'\ G is only approximate. 

In the A phase the effect is to require that the vectors d and 1 be parallel 
or antiparallel, and in fact by selecting one of the two configurations related 
by inversion, we can ensure that d = 1. In this case, we find 

H' A = 17(1), M' A = G'/H' A = SO(3). (56) 

In the B phase, the restriction is that R in (53) is no longer an unconstrained 
orthogonal matrix, but a rotation matrix through a definite angle (the 
Leggett angle 6*l = arccos(— \)) about an arbitrary axis n. Thus we find 



H' B = S0(2), M' B = G'/H' B = S 2 xS 1 . 



(57) 
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2.9. THE STANDARD MODEL OF PARTICLE PHYSICS 

There are remarkable similarities between the symmetry breaking pattern 
of 3 He and that found in the standard model of particle physics which in- 
corporates quantum chromo dynamics together with the unified electroweak 
theory of Weinberg and Salam. It is based on the symmetry group 

G = SU(3) col x SU(2) I x U(1) Y , (58) 

where I and Y denote respectively the weak isospin and weak hypercharge. 
The symmetry breaking from this down to the observed low energy sym- 
metry is described by the Higgs field, <fi, which plays the role of the order 
parameter. It is a two-component complex scalar field invariant under the 
colour group SU(3) co \, belonging to the fundamental 2-dimensional repre- 
sentation of SU (2)i, and with non-zero weak hypercharge Y = 1. It acquires 
a vacuum expectation value of the form 

(4>) = 4> = (° v ) , (59) 

thus reducing the symmetry to the subgroup 

H = SU(3) col x f/(l) em . (60) 

The generator of the remaining U(l) symmetry is the electromagnetic 
charge 

Q = h + \Y. (61) 

There may also be other stages of symmetry breaking at higher energies. 
The three independent coupling constants 53,52,51 corresponding to the 
three factors in G have a weak logarithmic energy dependence and appear to 
come to approximately the same value at an energy scale of about 10 15 GeV, 
especially if supersymmetry is incorporated into the model [11, 12, 13]. 
This suggests that there may be a grand unified theory (GUT) uniting the 
strong, weak and electromagnetic interactions in a single theory with a 
symmetry group such as SO(W). There would then be a phase transition 
(or a sequence of phase transitions) at that energy scale at which the GUT 
symmetry breaks to the symmetry group (58) of the standard model. If 
the model is super symmetric, then there must also be a supersymmetry- 
breaking transition. 

3. Defect formation 

The appearance of topological defects is a common feature of symmetry- 
breaking phase transitions. In this lecture, I shall review the defects asso- 
ciated with the various transitions discussed earlier, and the general condi- 
tions for the existence of defects. 
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3.1. DISCRETE SYMMETRY BREAKING 

The simplest possible field-theoretic model that exhibits symmetry break- 
ing is a model of a real scalar field described by the action integral 

1 = Jd 4 x [\{d^){d^) - ±A(0 2 - i? 2 ) 2 ] . (62) 

Here the action is invariant under the reflection symmetry cj) i— > — <f>. Thus 
the symmetry group is G = Z2, and the manifold of degenerate vacua 
reduces to a pair of points; the two vacuum states are characterized by 

<P = (4>) = ±V- (63) 

At high temperature, the equilibrium state is symmetric, with (</)) = 0. 
When the system cools through the critical temperature, (f> acquires a non- 
zero expectation value, but the sign is chosen arbitrarily. So it may happen 
that in one region, it chooses <fi ~ f] and in another <fi ~ —r]. When such 
regions meet, they must be separated by a planar defect, a domain wall, 
across which <p goes smoothly from one value to the other. The minimum 
energy configuration is determined by a balance between gradient energy 
and potential energy. At zero temperature one finds for example that a 
domain wall in the xy-plane is described by 

z 2 
(f>(z) = rytanh -, with £ = . (64) 
£ VAr/ 

As the system cools below T c , energy is trapped in the domain wall. In 
a sense the defect is a region of trapped old high-temperature phase, with 
the characteristic energy density that it had at T c . The wall is topologically 
stable. It can move, as one domain grows at the expense of the other, but it 
cannot simply break. A closed wall bounding a finite domain may of course 
shrink and eventually disappear. But this is a relatively slow process, so 
walls may have continuing effects. 



3.2. ABELIAN VORTICES OR STRINGS 

Now let us consider the case of an Abelian U(l) symmetry, such as that of 
superfluid 4 He. When the system is cooled through the transition temper- 
ature, the order parameter acquires a non-zero expectation value <p = ne 1 ^. 
The magnitude r\ is determined by the minimization of the free energy, but 
the phase d is arbitrary. It is chosen randomly. However, in a large system 
there is no reason why the same choice should be made everywhere; d may 
vary from one part of the system to another. We should expect the choice 
to be made independently in widely separated regions, especially if we are 
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talking about a transition in the early Universe, where such regions may 
have had no prior causal contact. 

When such a random choice is made, it may happen that around some 
large loop in space the value of -d varies through 2ir or a multiple thereof. 
In such a case, 4> must vanish somewhere inside the loop; indeed it must 
vanish all along a curve that threads through the loop. This is the core of 
a topological defect, a cosmic string or vortex. 

If the string is along the z axis, the order parameter around it typically 
takes the form 

cf>(r,<p,z)=r l f(r)e in v, (65) 

where r, tp, z are cylindrical polar coordinates, and n is an integer, the 
winding number. The function / has limiting values /(0) = 0, /(oo) = 1. It 
may be determined by minimizing the Ginzburg-Landau free energy, (27). 

For a superfluid, an important consequence of the expression (30) for 
the superfluid velocity is that the superfluid flow is irrotational: V A v s = 0. 
The vorticity vanishes everywhere, except in the core of the string, where 
the superfluid density vanishes. The string is a vortex. The form of (65) 
implies that there is a flow of superfluid around the string, with velocity 

(v 8 )<p = (66) 
m 4 r 

at large r, where m± is the mass of a 4 He atom. Thus the circulation around 
the string is quantized: 

<j> v s ■ dv = riK.4, (67) 
where the circulation quantum is 

k 4 = • 68 

There is a similar vortex in 3 He--B, but in that case is replaced by 
the mass of a Cooper pair, namely 2m^, so the circulation quantum is 

irh , . 

k 3 = — • 69 
m 3 

An important feature of the string or vortex is its topological stability 
deriving from this quantization. It can move around, but cannot break. 
A vortex loop can disappear by shrinking to a point, but a long, straight 
vortex is stable. 
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3.3. VORTICES IN A GAUGE THEORY 

If the symmetry is a gauge symmetry, with coupling to a gauge field A^, 
then around the string A has an azimuthal component, 

A ip (r,<p,z) = —g(r), (70) 
er 

where e is the charge and g has the same limiting values as /. By taking 
an integral round a large loop surrounding the string we find that it carries 
a quantized magnetic flux, 



/•27T 2irh 

$ = lim / A^rdip = n . (71) 

The magnetic field is given by 

B z = ^'(r). (72) 

(In the case of a superconductor e in the above should be replaced by the 
charge 2e of a Cooper pair, so the flux quantum is actually irh/e.) 

The functions / and g are determined by minimizing the free energy. 
For the Abelian Higgs model at zero temperature, with action integral (35) 
(setting c = h = 1), they satisfy the equations 



f" + - r f'-^(i-g) 2 f + W(i-f 2 )f = o, 
l 



g "__ g ' + 2 eWf{l-9) = 0. (73) 
r 

No analytic solution is known, but it is easy to find solutions numerically 
Note that there are two length scales governing the large-r behaviour of the 
functions, the inverse masses of the scalar and vector excitations (Higgs and 
gauge particles), 

m s 2 = 2Ar/ 2 , and m v 2 = 2e 2 ?? 2 . (74) 
The asymptotic behaviour depends on the ratio of these two, 

For k < 4 one finds that at large r 

1 - g oc r V2 e - mvr , 1 - / oc r- 1 / 2 e _m,r . (76) 
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In this case, the string has a narrow core that constitutes a magnetic flux 
tube, while the order parameter reaches its vacuum value over a larger 
distance. On the other hand, when k > 4, m v controls the behaviour of 
both 1 — / and 1 — g, with 

1 _ g oc r V2 e - mvr , 1 - / oc r- x e- 2m " r . (77) 

In superconductors [4], the two length scales are known as the correla- 
tion length £ = h/m s c and the Landau penetration depth A = h/m v c. Here 
large and small values of k distinguish so-called type-I from type-II super- 
conductors. In a type-II superconductor, vortices with |n| > 1 are unstable; 
there is a repulsive force between parallel n = 1 vortices which can stabilize 
a lattice of vortices. Hence there is an intermediate range of magnetic field 
strength within which the field penetrates the superconductor but confined 
to a lattice of flux tubes. 

3.4. DEFECTS IN NEMATICS 

It is easy to construct a stable linear defect similar to a string in the case 
of a nematic, by allowing the director n to rotate as one moves around the 
string through an angle tt: e.g., we can take 

n(r,(p,z) = (cos|,sin|,0), (78) 

provided we include in the order parameter Q a factor f(r) that vanishes 
at r = 0. Q then has no singularity because of the identification of n and 
— n. This configuration is called a disclination [7]. 

Like the superfluid vortex, the disclination is topologically stable. If the 
rotation angle were 2ir instead of tt it would not be. It could 'escape into 
the third dimension': at small r, we could rotate n upwards until at r = 
it points in the z direction, thus allowing /(0) to be nonvanishing. 

In addition to this linear defect there can be a point defect in a nematic 
liquid. Away from the centre, r = say, we can take 

n(r) = I, (79) 

again provided there is a factor in Q that vanishes at the centre. 

This is often called the hedgehog or monopole configuration. Like the 
vortex it is topologically stable, and cannot disappear spontaneously - 
though it may annihilate with an anti-hedgehog. 

3.5. THE FUNDAMENTAL GROUP 

The general conditions for the existence of defects can be expressed in terms 
of the topology of the vacuum manifold M, specifically its homotopy groups 
[14]. 
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The existence of linear defects for example is related to the possibility of 
finding non-trivial closed loops in M. Let us recall that different points in 
M correspond to different values of the order parameter labelling different 
degenerate vacua. A closed loop is a continuous map <fi : I — ► M from the 
unit interval of real numbers, I = [0, 1] C R, to M such that 0(0) = 0(1) 
(or equivalently a map from the circle S 1 to M) . Linear defects can exist if 
it is possible to find a closed loop that cannot be continuously shrunk to a 
point without leaving M, because then if the value of the order parameter 
around a loop in space follow this curve, it is not possible to fill in the 
values inside the loop continuously while remaining on M. 

In general, two closed loops are homotopic if it is possible to deform 
one continuously into the other within M. This is an equivalence relation, 
so we may define homotopy classes of loops. For example, when M is the 
circle S 1 , the homotopy classes may be labelled by the winding number, the 
(algebraic) number of times that we traverse the circle while going from 
to 1 in I. 

The homotopy classes constitute the elements of a group, the fundamen- 
tal group or first homotopy group of M, denoted by tti(M). To construct it, 
we introduce a base point b G M, and consider loops starting and finishing 
at b, i.e. 0(0) = (f)(1) = b. Then the product of two loops an d ip is the 
loop constructed by following and then ip: 



It is easy to show that this is a relation between homotopy classes, and so 
defines a product on the set of classes, which thus becomes the group n\ (M). 
The condition for the existence of topologically stable linear defects, strings 
or vortices, is that the fundamental group be non-trivial: tt\(M) / 1. 

For the Abelian case, where M = S 1 , the fundamental group is simply 
the group of integers, 7Ti(S' 1 ) = Z. The distinct possible linear defects are 
labelled by the elements of this group, the winding numbers. 

In a continuum version of the Heisenberg ferromagnet, we have M = S 2 . 
On the sphere, all loops can be shrunk to a point, so iri(S 2 ) = 1; there are 
no possible linear defects. 

For the nematic, however, the situation is different. On the sphere S 2 all 
loops are homotopically trivial, but this is no longer true when we identify 
opposite points to form RP 2 , because a curve that starts at one pole and 
ends at the opposite pole is closed in RP 2 but cannot be shrunk to a point. 
There is only one non-trivial homotopy class, because traversing the same 
loop twice gives a trivial loop; as we noted earlier, a disclination with a 
winding of 2tt rather than it is not stable. Hence the fundamental group in 
this case is tt\(RP 2 ) = Z2 = {0, 1}, the group of integers modulo 2. 




(80) 
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3.6. THE SECOND HOMOTOPY GROUP 

The conditions for the existence of other types of defects can also be ex- 
pressed in terms of homotopy groups. For point defects such as the hedgehog 
the relevant question is whether there are non-shrinkable two-surfaces in 



The second homotopy group [14] is defined in terms of closed two- 
surfaces, i.e. maps <\> : I 2 — > M from the unit square to M, such that 

0(0,*) = 0(1,*) =6, 0(s,O) = (f)(s,l) = b for alls and t. (81) 

In other words, 4> maps the whole boundary of I 2 to the chosen base point 
b e M. In effect, it is a map from S 2 to M, in which one designated point 
is mapped to b. 

Two closed surfaces are homotopic if one can be smoothly deformed into 
the other. This defines an equivalence relation, and hence a classification 
into homotopy classes of surfaces. 

As before, we can introduce a product on the set of closed surfaces, by 
setting 



This defines a product on the set of homotopy classes. (We could equally 
well have defined the product with the roles of s and t interchanged; it is 
easy to show that the results are homotopic to each other.) Thus we have 
defined the second homotopy group 7r2(M). 

As a simple example, let us consider a continuum version of the Heisen- 
berg ferromagnet, with SU (2) symmetry and an order parameter M trans- 
forming according to the 3-dimensional vector representation. Below the 
transition, the magnitude of M is fixed but its direction is arbitrary. Thus 
the subgroup H that leaves M invariant is H = £7(1) and M = SU(2)/U(l) 
= S 2 . In this case, the homotopy classes are labelled by an integer, the (al- 
gebraic) number of times the map wraps around the sphere. For example, 
a typical element of the homotopy class labelled by n is the map of one 
sphere on another defined in terms of polar coordinates 9, if by 



For this case, therefore, TT2(S 2 ) = Z, the group of integers. 

For the nematic, we have to identify opposite points of S 2 and pass to 
RP 2 but this makes no difference to the classification of closed surfaces. 
We again have TT2(RP 2 ) = Z, so the possible hedgehogs are labelled by an 
integer. 



M. 




(82) 



4> : S 2 -► S 2 : (0,<p) ' ^ (9, nip). 



(83) 
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One thing this classification cannot tell us, however, is whether con- 
figurations with |n| > 1 are actually stable. In some cases, it may be en- 
ergetically favourable for a configuration with winding number n = 2 for 
instance, to break up into two separate n = 1 configurations. Whether this 
actually happens is a question of detailed dynamics. 

3.7. DOMAIN WALLS 

As we noted earlier, domain walls occur when a discrete symmetry is bro- 
ken. More generally, the condition for the existence of domain walls is that 
the vacuum manifold M be disconnected. Domain walls are classified by 
the elements of what is often called the L zeroth homotopy group', denoted 
by ttq(M), whose elements are in one-to-one correspondence with the con- 
nected components of M. It is analogous to the higher homotopy groups: 
it may be regarded as classifying maps cf) : S° M, where S° is the 0- 
sphere (the boundary of the interval [—1,1] C R, namely the pair of points 
{1, —1}) in which the image of one chosen point is the base point of M. 

In a general case, the terminology is strictly speaking inaccurate, be- 
cause ttq(M) is not a group. There is one special case in which it is so, 
namely when H = 1, so that M is itself a group, M = G. In this case, 
the connected component Go of G containing the identity is an invariant 
subgroup (i.e., for any g G G, gG^g^ 1 = Go), and hence the quotient group 
G/Gq is defined; moreover 



Another way of characterizing ttq(M) is as a quotient of two groups. If 
(f> £ M can be connected by a continuous path to (fro, then one can always 
find a continuous path in G from the identity e to g such that D(g)<pQ = <f>. 
Hence the connected component Mq C M containing 0o may be identified 
with the set of elements {g<po : g £ Go}. The subgroup of Go which leaves 
0o unaltered is clearly H DGq. Hence, 



Now since Go is an invariant subgroup in G, it follows that H n Go is also 
an invariant subgroup of H. Thus H/{H n Go) is a group, and moreover a 
subgroup of G/Gq. One can then show that 



7T (G) = G/G . 



(84) 



M = G /(H nG ). 



(85) 



MM) = (G/G )/(H/(HnG )). 



(86) 



3.8. HELIUM-3 



Finally, let us return to the interesting case of 3 He, beginning with the 
superfluid 3 He--B phase, and initially ignoring the spin-orbit interaction. 
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We recall from (54) that M B = SO (3) x S 1 . Here both the zeroth and 
second homotopy groups tto(Mb) and ^(Mg) are trivial, so there are no 
topologically stable domain walls or monopoles. However, 

m{M B ) = Z 2 x Z, (87) 

so there are two different kinds of vortices. The factor Z classifies vortices 
around which the phase changes by 2mr, exactly as in the case of 4 He. 
However, the Z 2 factor arises because there are non-trivial loops in SO(3) = 
RP 3 . In a vortex corresponding to the non-trivial element of Z 2 there is 
no actual circulation around the string, but rather a relative rotation of 
the orbital and spin angular momenta. These are called spin vortices as 
opposed to mass vortices. 

Note that vortices may carry both types of quantum number simulta- 
neously. Such a combination is a spin-mass vortex. 

When we take account of the spin-orbit interaction the manifold is re- 
duced, according to (57), to M' B = S 2 x S 1 . In this case, we find 

Tr 2 (M' B ) = Z. (88) 

Viewed on a large scale there are monopole configurations. But since there 
are no short-range monopoles, these have no actual singularity. What hap- 
pens is that the order parameter near the monopole is forced to leave the 
manifold M' B , but can remain everywhere on the larger manifold Mb] the 
rotation angle in the order parameter, which is fixed to be the Leggett angle 
#L at large distances, can tend smoothly to zero at the centre, but (j) itself 
remains non-zero. 
We also find 

7n(Mk) = Z, (89) 

corresponding to the fact that the mass vortices are unaffected by the spin- 
orbit coupling, and survive to large distances. This is not the case, however, 
for the spin vortices, since there is no longer a Z 2 factor. What happens 
is that these become attached to a long-range soliton or domain-wall fea- 
ture. The order parameter at large distances around this vortex cannot lie 
everywhere on M' B , but to minimize the energy it does so except near one 
direction. Note that this is true in spite of the fact that ttq{M' b ) = 1, which 
means there are no truly stable domain walls: M' B is connected, but the 
relevant point is that it is not possible in Mb to deform the relevant loop 
in such a way that it lies entirely in M' B . (Such cases may be classified by 
the relative homotopy groups of Mb and its subspace M' B , in this case the 
group tti(Mb, M' b ).) 

Now let us turn to the A phase, for which according to (52) the 'vacuum 
manifold' is Ma = S 2 x S l O(3)/Z 2 . This space is again connected, so there 
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are no stable domain walls. However, we find 



7Ti(M A ) = Z 4 , 



7T 2 (M A ) = Z. 



(90) 



Thus there are monopoles, labelled by an integer winding number, and 
vortices labelled by an integer n modulo 4. On the other hand at long 
range the manifold, given by (56), is simply M' A = SO (3), whence 



Hence there are no stable monopoles and only one class of stable vortices; 
the latter are to be identified with the n = 2 short-range vortices. 

It is not hard to see what happens to the other short-range defects. 
For an n = ±1 short-range vortex, the corresponding loop in Ma cannot 
be deformed to lie entirely within M' A . In other words, we cannot make d 
parallel to 1 everywhere. The vortex becomes attached to a sheet or domain 
wall across which d rotates by it with a compensating rotation about 1. 

Similarly, around a short-range monopole we cannot deform the order 
parameter so that d remains everywhere parallel to 1. The monopole be- 
comes attached to a string in the centre of which d is in the opposite 
direction. 

4. Cosmology in the Laboratory 

Our present understanding of fundamental particle physics suggests that 
the Universe may have undergone a series of phase transitions very early 
in its history. One of the clearest signatures of these transitions would be 
the formation of stable topological defects with potentially significant cos- 
mological effects. To predict these we need to estimate how many defects 
would have been formed and how they would have evolved during the sub- 
sequent cosmic expansion. Calculations of the behaviour of the system in 
the highly non-equilibrium context of a rapid phase transition are problem- 
atic, however, and it is hard to know whether they are reliable. There is no 
direct way of testing them, because we cannot do experiments on the early 
Universe. 

But what we can do is to apply similar methods to analogous low- 
temperature transitions in condensed-matter systems, which often have 
a very similar mathematical description. Over the last few years, several 
experiments have been done in a variety of systems to test ideas drawn 
from cosmology. This has led to some extremely innovative and exciting 
condensed-matter physics, although the rsults are still somehwat confus- 
ing. 



*i{M' A ) = Z 2 



7T2(M' A ) = 1. 



(91) 
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4.1. DEFECT FORMATION IN THE EARLY UNIVERSE 

The electroweak transition, at about 100 GeV, where the W and Z particles 
acquire a mass through the Higgs mechanism, occurred when the age of the 
Universe was around 10~ 10 s. It is now believed that this is not in fact a 
genuine phase transition but rather a rapid but smooth crossover [15]. (This 
is possible only because this is a gauge theory.) There was probably a later 
transition, the quark-hadron transition at which the soup of quarks and 
gluons separated into individual hadrons. 

More interesting from a cosmological point of view, however, are the hy- 
pothetical transitions at even earlier times. If the idea of grand unification 
is correct, there would have been a phase transition of some kind at an en- 
ergy scale of around 10 15 GeV, corresponding to a time about 10 -36 s after 
the Big Bang. In some models, we expect a sequence of phase transitions, 
as the symmetry is broken in several stages, for example 

50(10) -► SU(5) -► SU(3) x SU(2) x U(l), (92) 

or 

50(10) -► SU(4) x SU{2) x SU{2) -► SU(3) x SU(2) x U(l). (93) 

The most attractive GUTs are supersymmetric. Since supersymmetry is not 
manifest at low energies, it must have been broken at some intermediate 
time, possibly yielding another phase transition, perhaps at about 1 TeV. 

Domain walls, strings and monopoles may all have been formed at early- 
Universe phase transtions, as indeed may more general composite objects 
of various kinds [8]. Monopoles and domain walls are cosmologically prob- 
lematic, for different reasons. Heavy domain walls, such as those that could 
have been formed in the early Universe, certainly do not exist in our Uni- 
verse today, and monopoles could be present only in very small numbers. 
So if these defects were produced at all, there must have been a mechanism 
to remove them completely or almost completely at an early stage. Inflation 
has often been invoked to do this job. 

Strings on the other hand could have survived in sufficient numbers to 
be cosmologically significant without violating any observational bounds. 
For a long time it was believed that they might serve to explain the ini- 
tial inhomogeneities in the density of the Universe from which galaxies and 
clusters later evolved. The idea that strings alone could seed these den- 
sity perturbations seems no longer viable, in the light of the data on the 
cosmic microwave background anisotropy. It is still perfectly possible to fit 
the data with models incorporating both strings and inflation [16], but a 
recent analysis concludes that strings probably do not make a significant 
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contribution [17]. They may, however, have had other important cosmolog- 
ical effects, for example in the generation of magnetic fields [18, 19, 20], 
high-energy cosmic rays [21, 22] and baryogenesis [23, 24]. 

For these reasons I shall restrict the discussion to the case of string 
formation. To be specific, let us consider the breaking of an Abelian £7(1) 
symmetry — though most of the discussion can easily be extended to non- 
Abelian symmetries. 

4.2. DEFECT FORMATION AT A FIRST-ORDER TRANSITION 

The nature of the early-Universe transitions is largely unknown, in partic- 
ular the order of each transition. In some cases, as I mentioned, there may 
be no true transition at all. Defects may be formed in any event, but the 
mechanism depends strongly on the order. In relation to condensed-matter 
analogues, most interest attaches to second-order transitions, and that is 
the case I will spend most time on. But I shall begin with what is in some 
ways the simpler situation of a first-order transition. 

In fact, the first 'cosmology in the laboratory' experiments were done 
with a first-order transition, namely the transition from normal isotropic 
liquid to nematic liquid crystal [25, 26]. 

Let us suppose, therefore, that there is a first-order transition, proceed- 
ing by bubble nucleation. Once the Universe reaches the relevant critical 
temperature, bubbles of the new low-temperature phase are born at ran- 
dom positions and start to grow until they eventually meet and merge. The 
nucleation rate 7 per unit space-time volume is given by an expression of 
the form 

7 (T) = A(T)e~ 5E(T) (94) 

where 5e is the Euclidean action for a tunnelling solution, and the prefactor 
A is typically of order T 4 . The nucleation rate determines the characteristic 
distance £ between nucleation sites, such that the number of separate bub- 
bles nucleating in a large volume V is U/£ 3 . Typically the bubbles expand 
at relativistic speeds, and then £ is of order 7 -1 / 4 . 

In each new bubble the order parameter <p becomes non-zero, and must 
choose a random phase There is no reason why there should be any 
correlation between the phases in different bubbles (except conceivably in 
the case of very near neighbours). So it is reasonable to assume that each is 
an independent random variable, uniformly distributed between and 2ir. 

When two bubbles meet, an equilibration process will occur, leading to 
a phase -d smoothly interpolating between #1 and #2 across the boundary. 
It is reasonable to assume that it will do so by the shortest possible path, so 
that the total variation will always be less than tt; this is called the geodesic 
rule. Numerical simulations have confirmed that it is usually though not 
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universally true — the rule may lead to a slight underestimate of the total 
number of defects formed [27, 28] . 

When these two bubbles encounter a third, it is possible that a string 
defect may be trapped along the line where they meet. This will happen 
if the net phase change from #i to -$2 to #3 and back to i?i is ±2tt rather 
than zero. If the geodesic rule applies and the three phases are strictly 
independent, the probablility of this happening can easily be seen to be |. 

Thus the total length of string formed in this process in a large volume 
V will be of order V/£ 2 . The length of string per unit volume will be 

L = ^, (95) 

where k is a numerical constant of order 1. (For example, if it is as- 
sumed that nucleation sites form a body-centred cubic lattice, one finds 
k = 3/2 7 / 6 = 1.34. A random lattice would really be more appropriate; 
that might well give a somewhat smaller value.) 

The first tests of this idea in condensed-matter systems were done in 
nematic liquid crystals, by studying the formation of disclination lines in 
the isotropic to nematic transition [25, 26]. The symmetry in that case 
is of course non-Abelian, but the principle is the same. We may assume 
that within each nucleating bubble of the nematic phase, the director n is 
an independent random variable, uniformly distributed over half the unit 
sphere (except near the walls where special effects come into play). 

The analogue of the geodesic rule is then the assumption that across 
the interface between two bubbles, the director always turns by an angle 
less than it. In that case, the probability that a disclination will be trapped 
between three bubbles with independently oriented directors is l/ir, so (95) 
should still hold. 

The experiments did in fact show reasonably good agreement with the 
predictions. Further experiments have been done to check the correlations 
between defects and antidefects [29]. 

4.3. SECOND-ORDER TRANSITIONS 

The argument is a little more complex in the case of a second-order phase 
transition. As the system cools through the critical temperature, the order 
parameter must acquire a non-zero value and choose a random phase. We 
may assume that the choice is made independently in widely separated 
regions. Thus there is a chance that defects will be trapped, and we should 
expect the formation of a random tangle of strings. What is less obvious is 
what the characteristic scale £ str of this tangle should be. Here £ str may be 
defined by the condition that in a randomly chosen volume £f tr there will 
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be on average a length £ str of string. In other words, the length of string 
per unit volume is 

£str = 75". (96) 
Sstr 

What determines £ str ? Obviously it is related to the correlation length 

of the order parameter, specifically of its phase. But this is not an answer. 
During a second-order phase transition, ^ is varying rapidly. Indeed, it is 
characteristic of second-order transitions that the equlibrium correlation 
length £ eq diverges at the critical temperature. So we must specify at what 
time or what temperature £ str should be compared with £ eq . 

An answer to this question has been given by Wojciech Zurek [30, 31, 
32], following an earlier suggestion of mine [33]. It is clear that in a real 
system going through the transition at a finite rate, the true correlation 
length £^ can never become infinite. In fact, for reasons of causality it can 
never increase faster than the speed of light. So, beyond the point where 
£eq = c, the adiabatic approximation, that ^ « £ eq (T), ceases to be valid, 
and instead one may assume that ^ will be more or less constant until 
after the transition, at least to the point where it again becomes equal to 
the decreasing £ cq . In a non-relativistic system, it is not the speed of light 
that is relevant, but some characteristic speed of the system. 

Zurek has given an alternative argument leading to essentially the same 
conclusion, based on a comparison of the quench rate and relaxation rate 
of the system. 

Let us assume that near the transition, the temperature varies linearly 
with time, so that 

T t 

e - 1 -f = f ( 97 ) 

J c ' q 

Here r q is the quench time. (We take t = when T = T c .) The equilibrium 
correlation length near T c has the form 

W r )=6H" v , (98) 

where v is a critical index. In mean field theory, v = ^, and this is often 
an adequate approximation. For 4 He, however, the renormalization group, 
gives a more accurate value, v = |. Similarly the relaxation time r diverges 
at T c : 

r(T) = T \e\-^ (99) 

where for 4 He the critical index \i = 1. This is the phenomenon of critical 
slowing down. The characteristic velocity is 

c(T) = = ^lel"-". (100) 

t{T) r 
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Note that it vanishes at T c . In 4 He, this is the speed of second sound, a 
thermal wave in which the normal and superfluid components oscillate in 
antiphase. 

Now information about the phase of the order parameter cannot prop- 
agate faster than the speed c(T). Hence after the transition the distance 
over which phase information can propagate is the sonic horizon 

Ht)= j; 4 T (mt - = YT l—^->. (101) 

This becomes equal to the equilibrium correlation length when 



ez 



(1 + M-^)- 



i 



(102) 



The time when this happens is the Zurek time 

*z = [(l + A*-")W] T k- ( 103 ) 

It is reasonable to suppose, at least as a first crude approximation, that 
at the Zurek time the characteristic length scale £ str of the tangle of strings 
or vortices should be equal to the correlation length: 

£str(*z) ~ & = Ccq(tz) ~ ?0 (^) • (104) 

Equivalently, we expect the density of strings or vortices (i.e., the length 
per unit volume) to be approximately l/£§, i.e., 



K / T 



2v 
1 + M 



W*z) = ^(jrJ , (105) 

where k is a numerical constant of order unity. Numerical simulations [34, 
35] suggest that it should in fact be somewhat less than unity, perhaps of 
order 0.1. Note that in 4 He, the exponent in (105) is \ in mean field theory, 
while using renormalization-group values it is | . This is the prediction that 
has to be tested. 



4.4. EXPERIMENTS IN HELIUM-4 

Zurek [30] initially suggested testing these predictions in superfluid 4 He. 
Experiments designed to test his predictions have been performed by Peter 
McClintock's group at Lancaster using a rapid pressure quench. The ex- 
perimental sample was contained in a small chamber that could be rapidly 
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expanded to lower the pressure, thereby sending it through the lambda 
transition into the superfluid phase. The number of vortices produced was 
found by measuring the attenuation of a second sound signal, generated by 
a small heater. 

The first experiment [36] did in fact see evidence of vorticity gener- 
ated during the quench, at roughly the predicted level. However, it was 
not conclusive for various reasons. Vorticity might have been produced by 
hydrodynamical effects at the walls. Also the capillary tube used to fill the 
chamber was closed at the outer end, so that during the expansion some 
helium was inevitably injected into the chamber, possibly again creating 
vorticity. Another problem was that it was not possible to measure the 
second-sound attenuation during the first 50 ms after the transition, so 
that later readings had to be extrapolated back to the relevant time. 

To overcome these problems, the apparatus was redesigned to minimize 
the hydrodynamic effects, and the experiment repeated [37]. Somewhat 
disappointingly, the result was null: no vorticity was detected with the im- 
proved apparatus. One possible explanation for this is that the vortices pro- 
duced may simply disappear too fast to be seen [38, 39]. The rate at which 
vorticity dissipates was measured in the rather different circumstances of 
vorticity generated by turbulent flow. It is not certain that the results can 
be carried over to the circumstances of this experiment. 

A third version of the experiment, incorporating further improvements, 
is now being planned [40]. Results are eagerly awaited. 

4.5. EXPERIMENTS IN HELIUM-3 

There are a number of advantages in using 3 He rather than 4 He. One is that 
because the correlation length is much longer (40 to 100 nm, rather than 
less than 1 nm), a continuum (Ginzburg-Landau) description is much more 
accurate than in 4 He. Moreover, the energy needed to generate a vortex is 
larger relative to the thermal energy, so it is easier to avoid extrinsic vortex 
formation. Another advantage is that since the nuclear spin is non-zero, 
one can use nuclear magnetic resonance to count the vortices. 

Perhaps the greatest advantage, however, lies in the fact that one can 
induce a temperature- rather than pressure-driven transition. This is be- 
cause of another characteristic of 3 He, namely that it is a very efficient 
neutron absorber, via the reaction 

n + 3 He^p + 3 H + 764keV. (106) 

Two experiments have been done with 3 He, one in Grenoble [41] and 
one in Helsinki [42]. Both use 3 He in the superfluid B phase, and look for 
evidence of vortices similar to those in 4 He. Both make use of the neutron 
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absorption reaction, by exposing the helium container to neutrons from 
a radioactive source. Each neutron absorbed releases 764 keV of energy, 
initially in the form of kinetic energy of the proton and triton. This serves to 
heat up a small region to above the transition temperature. It then rapidly 
cools, in a time of the order of 1 [is, and goes back through the transition 
into the superfluid phase. During this process we expect a random tangle 
of vortices to be generated. 

In other respects the experiments are very different. The Grenoble ex- 
periment [41], using a sample of 3 He-i? at a temperature much less than 
T c was essentially calorimetry. The total energy released, in the form of 
quasiparticles, following each neutron- absorption event was measured. Of 
the available 764 keV of energy about 50 keV is released in the form of ul- 
traviolet radiation. However, the measured energy was in the range 600 to 
650 keV, depending on the pressure, leaving a considerable shortfall. This 
is interpreted as being the energy lost to vortex formation. It is very hard 
to think of any other possible interpretation. 

The main feature of the Helsinki experiment [42], using a sample of 3 He- 
B at a considerably higher temperature, not far below T c , was the use of 
a rotating cryostat. If a container of helium is rotated rapidly, vortices are 
generated at the walls and migrate to form a central cluster parallel to the 
rotation axis. However, if the rotation is slower, no vortices can be formed. 
In 3 He--B, it is possible to ensure that no vortices at all are present. We 
then have a remarkable situation. The normal fluid component is rotating 
with the container, but the superfluid component, which cannot support 
vorticity, is completely stationary. Thus there is a counterflow velocity, a 
difference v = v s — v n between the velocities of the two components. This 
introduces novel hydrodynamic effects; in particular a superfluid vortex 
moving relative to the normal fluid experiences the transverse Magnus force. 

In consequence vortices above a certain minimum size tq and correctly 
oriented are expanded until they reach the walls of the container, and then 
migrate to join a central cluster parallel to the axis. The number of vortices 
'captured' in this way can be determined by nuclear magnetic resonance 
(NMR) measurements. It is possible to detect each individual vortex joining 
the cluster. 

The number of vortices we expect to be captured can be predicted. It is 
essentially the number of vortices with sizes between the required minimum 
size ro, which depends on the counterflow velocity v, and the maximum ra- 
dius of the bubble. The size distribution of loops formed is expected to be 
scale invariant. This leads to a very simple prediction. There is a critical ve- 
locity v cn for neutron-induced vortex formation, which is substantially lower 
than the critical velocity v c for spontaneous vortex formation at the walls. 
If v > v^, the number of vortices captured after each neutron-absorption 
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event should have the form 




(107) 



where c is a calculable constant. Remarkably enough, all the dependence on 
the bulk temperature, the pressure and the magnetic field is contained in 
the value of v cn . Hence if N is plotted against v 3 for various values of these 
parameters, one should see a set of straight lines with a common intercept 
at — c on the vertical axis. This simple prediction does fit the experimen- 
tal results very well over a considerable parameter range, providing good 
evidence for the validity of the prediction. 

It has also been possible to test the predicted dependence of v cn on 
temperature, namely v cn oc e 1 / 3 . This again is a good fit to the data. 



4.6. EXPERIMENTS IN SUPERCONDUCTORS 

It is particularly interesting to test the predictions of defect formation in 
superconductors, because they provide an example of a gauge theory. 

The first experiments [43] were done by a group at Technion, using a thin 
film of the high-temperature superconductor YBCO. The film was raised 
above the critical temperature by shining a light on it, and then allowed to 
cool. The object of the experiment was to determine the number of defects 
formed, in this case 'fluxons' each carrying one quantum of magnetic flux. 
What Carmi and Polturak measured, using a SQUID detector, was actually 
the net flux, i.e., the difference AiV = N + — iV_ between the numbers 
of fluxons and antifluxons. In fact they saw no evidence for any fluxon 
formation, with an upper limit of |AiV| < 10. 

This result has to be compared with predictions based on Zurek's work. 
In this case the Zurek length £z is estimated to be about 10 -7 m, so within 
the 1 cm 2 sample we should expect the total number of defects to be 

N = N + + iV_ « 10 10 . (108) 

The net flux may be estimated by assuming that the phase of the order 
parameter performs a random walk with a step length of £z and a typical 
angle 5 ~ it/ 2. This suggests that 

where L ~ 20 mm is the perimeter of the sample. (Note that according to 
this argument AiV is of order iV 1 / 4 .) This yields 



AiV « 100, 



(110) 
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in clear contradiction to the experimental results. It should be noted that 
this prediction is based on (105) with the constant k set equal to unity, so 
there may be scope for reducing it slightly, though probably not by enough 
to remove the discrepancy. 

On the other hand, Carmi and Polturak in fact suggest that the dis- 
agreement is more serious, because in a gauge theory the mechanism of 
defect formation is different and the geodesic rule is unreliable [44, 45, 46], 
so one should perhaps expect AN to be of order N 1 / 2 , leading to an esti- 
mate AiV k, 10 4 which is obviously in very severe disagreement with the 
results. This is a point that needs further theoretical study. 

However, the same group have also performed another experiment [47], 
with very different results. This involved a loop of superconducting wire laid 
down in a square-wave pattern across a grain boundary in the substrate so 
as to create a series of N = 214 Josephson junctions in series. As the wire 
cools it becomes superconducting before the Josephson junctions start to 
conduct, so in effect each segment of wire between neighbouring junctions 
is initially a separate system, so it is reasonable to assume that their phases 
are random and uncorrelated. Hence some flux will be trapped when the 
wire eventually becomes a single superconducting loop. The experiment 
revealed an r.m.s. flux of 

AiV cxp = 7.4 ± 0.7. (Ill) 
The theoretical prediction in this case would be 

AiV th = ivW = 3.6. (112) 

It is perhaps rather surprising that the experiment saw more flux than 
predicted. The authors suggest that this may again be due to a breakdown 
of the geodesic rule with an r.m.s. value of 5 closer to tt than to ir/2. 
(Arguably, if 5 is uniformly distributed between —tt and tt, we should use 
an r.m.s. value of tt/-\/3 rather than 7r/2, but the difference is minimal, 
leading to AAT th = 4.3.) There could also perhaps be a non-zero phase 
change along the section of the loop away from the Josephson junctions. 

Recently experiments have been performed by a different group [48, 49] 
on annular Josephson tunnelling junctions, comprising two rings of super- 
conducting material separated by a thin layer. When the system is cooled 
through the critical temperature and the rings become superconducting, 
one may expect that the random choice of phase will lead to trapping of 
fluxons. For the experiments done so far the predicted number trapped is 
less than one fluxon on average, which is not ideal. Nevertheless, they have 
detected flux trapping at roughly the predicted level. An important feature 
of this experiment is that it is possible to vary the quench rate and so test 
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the dependence of the number of fluxons on the quench rate r q , as given 
by (105). The results are consistent with the Zurek predictions, though the 
scatter is large. 

4.7. DISCUSSION 

Experiments with 3 He, with liquid crystals and with superconducting loops 
have all confirmed the basic idea that defects are formed during rapid phase 
transitions. The best evidence so far that Zurek's predictions of defect num- 
bers are sound comes from the 3 He experiments, though the others are 
reasonably consistent. 

On the other hand, neither the 4 He experiment nor that with a super- 
conducting film have shown any evidence for defect formation. 

At first sight, the discrepancy between the results with 4 He and 3 He 
may be surprising, but in fact the differences between the two systems are 
very great. Karra and Rivers [50] have argued that a very important factor 
is the great discrepancy between the widths of the 'critical region', below 
the critical temperature and above the Ginzburg temperature Tq [51]. This 
is the temperature above which thermal fluctuations are large enough to 
create a significant transient population of thermally excited small vortex 
loops. It is given approximately by the condition that 

£ 3 q (T G )AF(T G ) = k B T G , (113) 

where AF is the difference in free energy between the 'false- vacuum' state 
with = and the broken-symmetry equilibrium state. Above this temper- 
ature, it appears, the formation of long vortices is suppressed. It happens 
that 3 He and 4 He are very different in regard to the width of the critical 
region between Tq and T c . In 3 He it is extremely narrow; Tq is very close to 
T c , at e ~ 10~ 8 . In 4 He, on the other hand, Tq is about half a degree below 
T c . Karra and Rivers [50] used thermal field theory, with a Gaussian ap- 
proximation, to show that the Zurek predictions should be approximately 
valid provided that 

<T G )^ < 100, (114) 

TO 

a condition that is very well satisfied for the 3 He experiments where the 
left hand side is about 10~ 5 and badly violated for those in 4 He, where it 
is 10 10 . 

Also puzzling is the discrepancy between the different experiments in 
superconductors. There is some doubt about how to compute the number 
of defects formed in a transition in a theory with a local gauge symme- 
try. There is another mechanism operating in a gauge theory [46, 52], but 
if anything this makes the discrepancy more puzzling because it tends to 
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suggest that the Zurek prediction of defect numbers would be an underesti- 
mate. On the other hand, it is worth noting that the inequality (114) is also 
seriously violated in the superconducting film experiment, though whether 
the argument leading to it is valid in the case of symmetry breaking in a 
gauge theory is not clear. 

What is clear is that there is as yet no certainty about when the 
cosmology-based predictions of defect numbers are reliable. Only further 
experimental and theoretical work will resolve this question. 
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